LAPLACE INVARIANTS FOR GENERAL HYPERBOLIC SYSTEMS 



CHRIS ATHORNE AND HALIS YILMAZ 



Abstract. We consider the generalization of Laplace invariants to linear differential sys- 
tems of arbitrary rank and dimension. We discuss completeness of certain subsets of invari- 
ants. 



1. Introduction 

The classical Laplace invariants j2] were introduced in the context of second order, linear 
hyperbolic systems of the form 

z, xy +az, x +bz, y +cz = (1.1) 

where a, b and c are given functions and z = z(x, y) is an unspecified solution of this partial 
differential equation. 

The form of equation (jl.ip is unchanged under a general transformation z h-> g(x,y)z 
where g{x,y) is a sufficiently differentiable, but otherwise arbitrary, function. In fact the 
coefficients of the equation are simply mapped into new functions, 

a a' = a + g~ x g, y , 
b h-> b' = b + g~ x g, x , 

c d = c + g~ 1 ag, x +g~ 1 bg, y +g~ 1 g, xy (1.2) 

and it is easily seen that the following two functions are invariant under such a transforma- 
tion: 

h = a, x +ab — c, (1.3) 
k = b, y +ab-c. (1.4) 

More than this, the pair {h, k} is a complete set of invariants in that two equations of 
the form f 1 1.11) having exactly the same invariants, as functions of x and y, must necessarily 
be related by a gauge transformation of the sort described. The family of equations is thus 
partitioned into equivalence classes labelled by these pairs of functions. These functions are 
called Laplace invariants by many researchers in integrability theory (see e.g. [3], [I], [8], [9], 

TO- 

Such invariants have played an important role in recent work on the geometrical theory of 
integrable systems and soliton equations. It is not our purpose to rehearse these connections 
here and we refer the interested reader to references [HI [32] where much of the material 
is reviewed. However, it is important to point out that a valuable role is played by the 
Laplace map, a differential map between equations of the form (II. ip which acts on the 
equivalence classes according to the equations of the two-dimensional Toda lattice [TT] . 
The generalization of the Laplace map to higher dimension and higher rank systems is of 
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great importance [TJ [12]. This paper should be regarded as a prologomenon to a general 
theory of such transformations. 

Before proceeding let us note that the form (11. ip . though symmetric, has a degree of 
redundancy about it. We may choose to transform it using a gauge transformation z h-> gz 
where g satisfies g, y = —a(x, y)g. In this case the transformed equation is 

z, xy + j {k - h}dy.z, y -hz = 0, (1.5) 

and the dependence on the equivalence class is explicit. An equation of this reduced form, 

z, xy +bz, y +cz = 0, (1.6) 

still retains a gauge covariance, namely z i— > g(x)z, the gauge function depending upon x 
alone and it is naturally written as a system in z and z, y : 

d x —(3c \ I —(3z, y j n q 



(1.10) 



1//3 d y J \ z 

where (3, x = (3b. 

Of course, we might equally consider reduced forms 

z,xy +az, x +cz = 0, (1.8) 

with y dependent gauge transformations, but what we cannot do in general is reduce to the 
form 

z, xy +cz = 0, (1.9) 

as this requires that the special relationship h — k, Wx, y should hold. 
Equally we could start with a general system form 

d x + /in h n \ f zt 
h 2 \ d y + h 22 J V ^2 

as is done in pQ. Gauge transformations preserving this form of system are 2x2 diagonal 
matrices acting on the two component vector of the Z{. The gauge invariants are 

(12) = h l2 h 2U (1.11) 

[12] = h u ,y-h 22 , x +l\n(^), X y. (1.12) 

Z n 21 

However the redundancy is also present here and we can use the gauge transformation to 
kill the diagonal terms h\\ and h 22 . This leaves us with the canonical form 

dx h 12 \ f z x 
hi d y ) \z 2 

and residual gauge transformations 



0, (1.13) 



with invariants 



(12) = h 12 h 21 , (1.15) 

[12] = iH^U- (1.16) 
z n 2 i 
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It is not difficult to verify that these invariants are a complete set for the canonical form 

dm. 

In what follows we shall consider n x n systems and discuss the completeness of the sets of 
invariants constructed in a similar manner to those presented in this introduction. We shall 
also relate them to second order, matrix equations, i.e. those of the type ( II. ip but having a, 
b and c as square matrices rather than simple functions. 

We use the word dimension to denote the number of independent variables which we shall 
henceforth write as X\ , %2 ■•• ; X n . By rank we shall understand the number of components 
in the solution vector z: z\, z%, . . . , z r . 



h 



21 



d 2 + h 2 2 



hln 
h2n 



\ 



d + h I 



2. Invariants For General Hyperbolic Systems 
Definition 2.1. Let L be an n x n matrix differential operator 

/ di + /in hi2 

L = 

where di stands for d/dxi and the hij are functions of x±, X2,.--,x n . If g is a diagonal n x n 
matrix such that g~ x exists, then H = H(hij) is invariant under the gauge transformation 

L' = g- l hg, 

so long as H = H(hij). 

2.1. The case where rank and dimension are equal. In this case we deal with matrix 
differential operators 



L 



/ <9i 

do 



\ 

and gauge transformations 
of the form 



\ 
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h\2 
h 2 2 



h<2n 



hnn ) 



L !->• L' = g^JLg 



( 9i 

g 2 
\ ... 








9n J 



The and gi here are functions of all variables xi,X2, ■ ■ ■ 
reduced (canonical) form in which the diagonal entries hn, h 2 %, 
solving the n equations: +hugi — 0, i — 1, 2, . . . , n. 



(2.1) 



(2.2) 



(2.3) 



x n but we may choose the 
. . . , h nn are gauged away by 
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/ dx 

o d 2 
V o ... 
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(2.4) 



The residual gauge freedom is 
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( 9i(xi) 
\ 








(2.5) 



• • • 9ni&n) J 

where hatted variables are deleted from the list of arguments in each g^. Under such trans- 
formations 

hj ^ 9i{xi)~ l 9j{xj)hij (2.6) 
and it is easily seen that the following objects are all invariant: Choose from the n labels 
{1,2, ... ,n} a subset of p distinct ones, {i\, i 2 , ■ ■ ■ , i p }, and define the symbol: 



{1X12 ...l p ) 



hi x i 2 hi 2 i 3 . . . hi p i 1 . 



(2.7) 



We say the symbol {i\i2 ■ ■ ■ i P ) has length p. Thus in the case of the symbols of lengths 2 and 
3 we have (ij) = hijhji and (ijk) = hijh^h^. 

Because of the cyclic symmetry in these products there will be p ^- P y. sy m bols of length 
p. The symbols of length p are permuted under the action of S n , the symmetric group on n 
labels. 

In addition there are \n{n — 1) invariants denoted by square bracket symbols thus: 



(2.8) 



We call the invariants (12. 7ft . (12. 8ft simple. All functions of these symbols are themselves 
invariant but we will now show that within the set of simple invariants there are a complete 
subset i.e. a set the knowledge of which is enough to determine the operator L completely 
up to gauge transformations. 



Lemma 2.2. The functions [ij] and (i^ ■ ■ - i p ) are invariants. 
Proof. We consider the n x n differential operator matrix L 

• h ln \ 

■ 



L 



/ di hi2 
h 2 \ d 2 



\ h n i h 



n2 



n J 



where hij are functions of x\, X2,..., x n . We find the invariants of L by using the gauge 
transformation, g^Lg = L', where g is a n x n diagonal matrix 



(9i 

92 

V 



Then g x Lg = L' gives us 



(ln#i),i> 



\ 




9n J 



(i ^ h 



(2.9) 
(2.10) 



Now 

Z Ilji Z llji Cji 

which gives 
since 

ft., P = 0, (2.11) 
where r = This gives us the antisymmetric invariants 

M = \Q*T L )m- (2-12) 
Finally we consider the following relations 

hi 3 i4 = 9% z 9iJ l iz%^ 



Then we obtain 



hi p -ii p 9i p _ 1 9i p hi p _ 1 i p 

% p %\ *?tp O il lpi\ 



h'- ■ h'- ■ h'- 1/ . h'- ■ — h ■ ■ h ■ h ■ h h 



to give the p-index invariants: 

(^1^2^3 • • • l i"p) hi 1 i 2 hi 2 i 3 hi^i A . . . h{ p—iiphipiu (2-13) 

where the i r are a choice of p distinct integers in {1, 2, . . . , n}. 

By recalling f)2.12p and ( I2.13P we now collect all the invariants of L as follows: 

(^i^2^3 • • • ip) hi 1 i 2 hi 2 i 3 hi^i i . . . hi p _ 1 i p hi p i 1 
Definition 2.3. The functions [ij] and (i^is • • -i p ) are called the simple invariants of L. 
Theorem 2.4. The simple invariants form a complete set for the equivalence class of L 



under gauge transformations, where L is defined by (2.4)- 

Proof. The proof depends on showing that one can construct a suitable gauge matrix g. 
In other words we need to show that 

M = M 
(iii 2 i 3 ■ ■ ■ i P )' = 0M2«3 ---i P ) 

where {i h i 2 , £3, - - - , i p } C {1, 2, . . . , n}. 



U = g~ x hg 



We already know that '=!>■' is true. We only need to prove the part. Assume the RHS 
is true i.e. 

M' = M 
(iii 2 i 3 ■ ■ ■ ip)' = {hi 2 h •• -i P ), 

,i P } Q {1,2,..., n}. Let us choose amixn diagonal matrix / such 



for all subsets i 2 , 
that 



where 
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( h o ... o \ 
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V o o ... /J 
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■ K- 
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In 
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= h 12 h 23 h34 ■ 


■ K- 
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Then we obtain 



where 



fn-l — ^12^23^34 • • • h n -2n-lhn-ln } 
Jn — ''12 "-23 "-34 • • • u n-2n-X u n-Xn- 



/ dx + /in ^12 
h 2 i d 2 + h 22 



nl 



h 



n2 



Thus we need to show 



We easily prove f )2.14p as follows 



hln \ 

h 2n 
d + h J 



= (In/O,* (z = l,2,...,n). 



(2.14) 



Let i < j. Then 



/i 
f) 



^12^23 • 

h 12 h 23 . 



h' h' 



■ hj-ijhjj+i 



■ ■ fort— In 
• • fori— In 



Thus 



Ji 



h' h' h' 

' l ii+V L i+li+2 ■ ■ ■ "'j-lj 

hu+lhi+ii+2 ■ ■ ■ hj-ij 



h' h' 



tij-ijh'ji hji_ h 



hu+ihi+ii + 2 ■ ■ ■ hj-ijhji h'ji 

(ii + li + 2...j)' (ij) Kj_ 
(n + li + 2... 3 )- h' 3l -h> t3 

(ij)'' « 



l 3 



K 



since (ii + li + 2... j)' — (ii + li + 2 
Similarly 



hji = ^-h 



.j) and (ij)' = (ij). 
hu+ihi+ii+2 ■ ■ ■ hj-ij 

= y y y ji 

' b ii+l' b i+U+2 ■ ■ ■ ,L j-lj 

(n + li + 2... j) 
(ii • 1/ • •_'..../)' " 



= tin- 



Hence for i 7^ j we obtain 



So we have 



( d 1 + h u 



h' 



hi 
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12 _ 
d 2 + h 2 2 
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~l~ h nn J 



where 

hu = (m 

We now need to seek a single function 9 so that 

9-\f- 1 Lf)9 = U. 
This requires that 9 satisfy the following equations: 

9~ 1 9 ii + h u = 

i.e. 



9 : i — —hu9. 



The above equations are consistent 



,v ,3 



(9 tj ) ti , which gives 



h 



h 



(2.15) 



'3 3, i 



(2.16) 



(Si 
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hji 


fi 






fi 
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Recalling (I2.15P we write 

ha = (ln/i) ti , 

and if we substitute these into the equation (12.161) we obtain 

M' = M 

since 



where 



So the equality of invariants guarantees that the Frobenius integrability condition is satisfied: 
there exists a function 9 such that 6~ l (f~ ll Lf)6 = L' i.e. 

g- l hg = L\ 

where g = Of . Hence the given invariants of L are a complete set. 

It should be noted that the simple invariants are not algebraically independent. For 
instance, 

(ijk)(ikj) = (ij)(jk)(ki) (2.17) 
so that there must be a smallest set of simple invariants which is still complete. A minimal 
complete set is given in the following result: 

Theorem 2.5. The simple invariants (lz) ; [ij] and (lij) form a minimal complete set. 

First we prove some lemmas. 
Lemma 2.6. We consider a simple invariant of length m 

(il?2^3 • • • im—l^m) ^ill2 ^«2i3 ^«3«4 ■ ■ ■ ^'im-lim^'imil ■ (2.18) 

Let m be a positive integer such that m > 4. Then 

f \ («l*2*3---«m-l)(*l*m-l«m) , n ln x 

(lll2*3 • • • «m-l«m) = J—- v (2-19) 

(Vm-l) 

Proof. 

DCJO {kkk- ■ ■im-x){Hi m -xim) 
ixri o = r 

(«iVi) 

_ hjrhJ^WS ■ ■ ■ him-2irn-\hi m _ 1 i 1 -hi 1 i m _ 1 hi rn _ 1 i m hi rn j il 

hi 1 i m _ 1 hi m _ 1 i 1 
— hi-Lizhfoiz . . . hi m _ 2 i m _ 1 hi m _ 1 i m hi m i 1 
= (iii 2 i3 ■ ■ ■ im-iim) = LHS 

Hence we can replace simple invariants of length m > 4 with invariants of length m — 1 up 
to multiples of invariants of lengths 2 and 3. 



Lemma 2.7. Let i,j, k be three positive integers such that i^j^k . Then 



Proof. 



Similarly 



(V) 
(ijk) 

(iO(ij) 

(iij)(ijk)(iki) 
(i<)(ij)(i*) 



(iO(ij) 

(iij)(ijfc)(ifci) 

(lz)(lj)(lfc) 

huhijhji.hijhjihii 

huhix-hijhji 
hijhji = (ij) = LHS 

hiihijhji.hijhjkhki.hikhkiha 
hiihu.hijhji.hikhki 

hijjjkhki 

(ijk) = LHS 



(2.20) 
(2.21) 



Lemma 2.8. The invariants (lij) are irreducible (i.e. they cannot be written purely in terms 
of invariants with length 2) 

Proof. We will prove this by contradiction. So assume (lij) is reducible. Thus (lij) can 
be expressed in terms of the invariants (lz), (lj) and (ij). So let 

(Uj) = F[(U),(lj),(tj)} (2.22) 

If we differentiate the equation (I2.22p with respect to ha, h\j and hji respectively we obtain 
the following partial differential equations: 









d(lij) 


dF 


dha 


d(li) 


d(lij) 


dF 


dh^ 


d(lj) 


d(lij) 


OF 


dhji 


d(ij) 



since (lij) = huh^hji is independent of ha, h\j and hp 
Thus, we find 

dF 



d(li) 
dF 

Wf) 
dF 
d(ij) 









since hu ^ 0, hji ^ and ^ 0. 



This shows that (lij) = constant. This is a contradiction. Therefore the invariant (lij) 
is irreducible. 

Proof of Theorem \2.5l We have considered the following simple invariants of length m: 

First we have shown (Lemma 12.61) that these invariants can be reduced up to length 3 and 
then we have shown (Lemma \2.7\i that the invariant (ij) can be written in terms of the 
simple invariants (li) and (lij) and we have also proved that the simple invariant (ijk) can 
be expressed in terms of the invariants (li) and (lij). Finally, we have proved (Lemma 
12. 8p that the invariant (lij) is not reducible, in other words, it can not be reduced to the 
invariant of length 2. 

Hence the proof of the theorem is complete and the result follows: Any invariant of length 
m can be written in terms of the minimal invariants (li) and (lij) where these minimal 
invariants together with [ij] form a complete set. 



3. Matrix Covariants For General Hyperbolic Systems 
3.1. Matrix Covariants. Let us consider the system 

Lz = (d x d y + ad x + bd y + c)z = (3.1) 

where a, b and c are m x m square matrices. This case is considered in [5]. The gauge 
transformation on the differential operator L is V = g~ 1 Lg, where g is a m x m diagonal 
matrix which gives 

h = a T + ba — c 

, l l ( 3 - 2 ) 
k = b t y + ab — c 

where h and k are gauge covariants for the system fl3.ll) : h! = g _1 hg, k' = g~ l kg. These 
covariants are sometimes called invariants in the literature [5]. 



3.2. Matrix Covariants for L. Let us consider L as a (ni\ + ni2) x (m^ + mi) differential 
matrix operator such that 

\ h 2 i o 2 + h 22 J 

where h 1± e M mimi , h 12 G M mim2 , h 21 G M mami , h 22 G M m2m2 and M mimj is the set of 
rrii x rrij matrices. 

Strictly speaking we should write I mi d\ and I m2 d 2 for the differential operator entries 
where I mi , I m2 are the unit matrices of dimensions m x and m 2 . This should be understood 
in what follows. 

The 'gauge' transformation g on L is L' = g~ l hg for 
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where g\ G M mimi and g 2 €= M m2ni2 are both invertible square matrix functions of Xi,X2- 
Under this action, L' = g^hg, we have 

Ki = 9i 1 hugi + g^gi,!, 
K2 = 92 1 h2292 + g 2 ~ 1 g2,2- 

3.3. Definitions. We call an object H of type G, t x Gj if H' = g~ l Hgj. Therefore h\2 is of 
type Gi x G 2 and h 2 i is of type G 2 x Gi. Covariants are of type Gj x Gj. In other words if is 
a covariant if if' = g~ x Hgi. Invariants are given by the traces of covariants. The operators 
d\ + /in and d 2 + h 22 are of types G\ x Gi and G 2 x G 2 respectively: 

9i + /i'n = gi 1 {di + h ll )g l , 
d 2 + h' 22 = g 2 ~ 1 (d 2 + h 22 )g2- 

But they are differential operator covariants . We seek matrix covariants . The simplest 
matrix covariants are ^12^21 of type Gi x Gi and ^21^12 of type G2 x G2 since, 

h 21 h 'l2 = 92 1 ( h 2ihi2)g2- 

Let us call (12) = hi 2 h 2 i and (21) = Zi2i^i2, where (12) e M mimi and (21) e M m2rrt2 . Note 
that we use similar notation to before but now (12) ^ (21). Our aim is now to form higher 
matrix covariants. For simplicity we call d\ + h n = D\ and d 2 + ^22 = D 2 . The operators 
Di and D 2 are of type G x x Gi and G 2 x G 2 respectively. Therefore one easily see that 

h,2\D\ and -D2^2i are of type G 2 x Gi, 
hi 2 D 2 and are of type Gi x G 2 . 

Hence 

di = h 12 D 2 h2iDi - D 1 hi 2 D 2 h 2 i is type of Gi x Gi, 
C22 = h 2 iD 1 h 12 D 2 - D 2 h 2 iD l h 12 is type of G 2 x G 2 . 

But these are still not matrix covariants, since they have leading differential operator terms 

c 11 = -[A,(12)]^ + ... 

C22 = -[D 2 ,(21)]d 1 + ... 

We would like to subtract off multiples of <9 2 + /i 22 from en and d 1 + h u from c 2 2 to remove 
the differential operators but each operator is of the wrong type. To circumvent this we turn 
Cn, C22 into respectively G2 x G2 and Gi x Gi of type covariants by: 

h 2 iCuh 12 = -h 21 [D 1 , (12)] d 2 h 12 + .. .matrix 

= -h 21 [D u (12)}h 12 (d 2 + h 22 ) + ■ ■ ■ matrix, (3.3) 

h 12 c 22 h 2l = -h 12 [D 2 , (21)]h 21 (d 1 + h n ) + . . . matrix. (3.4) 
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Since each part in expression (I3.3P is now of type G 2 X G 2 and each part in (13. 4p of type 
G\ x G\, we must have matrix covariants: 

[12] = h 12 c 22 h 21 + h 12 [D 2 , {21)}h 21 {d 1 + h n ), 
[21] = h 21 c n h 12 + h 21 [D u (12)]h 12 (d 2 + h 22 ). 

Simplifying these give 

[12] = (12) [D u h 12 D 2 h 21 ] - h 12 D 2 h 21 [D u (12)] of type G x x G 1 (3.5) 
[21] = (21)[D 2 , h 2l D x h 12 ) - h 2X D x h X2 \D 2 , (21)] of type G 2 x G 2 (3.6) 

as matrix covariants where [12] G M mimi and [21] G M mam2 . 
The case m x = m 2 = \: 

We find a reduction of [12] and [21] in the case m\ = m 2 = 1. So in this case (12) = h\ 2 h 2 \ 
and (21) = h 2 \h\ 2 are just equal functions and (21) = (12) = (12), the earlier invariant. By 
substituting D\ — d\ + h\\ and D 2 = d 2 + h 22 in the covariants [12], [21] and then by doing 
some differential and algebraic calculations we obtain the function covariants as follows: 

[12] = -^(12) 2 12 -(12) 2 [12], 

\ (3-7) 
[21] = - i (12) 2 12 + (12) 2 [12]. 

One easily sees that 

[12] + [21] = -I(12)f u , m 
[21] - [12] = 2(12) 2 [12] 

Thus relating the expressions from the new covariants to the old invariants in this case 
(rrii — m 2 — 1)- 

3.4. The case where rank exceeds dimension. It is clear that in the case where the 
rank r is larger than the dimension n we may attempt to repeat the arguments of section 
|2] under the weaker hypothesis that the and gi are matrices and no longer (commuting) 
functions. The canonical form (12. 4p still suffices where now the hij are rectangular matrices 
of type rrii x rrij, where an m 8 x m 8 unit matrix is taken to stand (but omitted) before each 
operator, <9j, and where m% + m 2 + . . . + m n = r. 
The case n = 2: 

In this case we consider a differential matrix operator L such that 

L=ff! h ?) (3-9) 



x hi d 2 / 

where h\ 2 G M mirn2 and h 2 i G M m2mi are matrix functions of x\ and x 2 . We have assumed 
a gauge transformation to this form as before. 
The gauge transformation 

L^L' = g^hg, (3.10) 

where 

9 = ( 9li * 2) ? (3.11) 
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9i 

92 



gives us 

hi 2 H> h' 12 = g^hvigi (3.12) 

h 2 \ i ^ h 2 i = 92 lh 2i9i (3.13) 

where #i and #2 are invertible square matrices such that g\ G M mimi and g 2 G M mam2 . 
So the relations (I3.12p and (13.131) give us 

= S'r 1 ( h i2h 2 i)gi 

K1K2 = 92 1 (hih 12 )g 2 

Thus we have 

(12) / = ^r 1 (12)^ 1 (3.14) 

(21)' = g 2 \2\)g 2 (3.15) 

where (12) G M mimi and (21) G M m2ni2 are matrix covariants such that 

(12) = h 12 h 21 (3.16) 

(21) = h 21 h 12 (3.17) 

By doing some algebraic calculations over (I3.12p and (13.131) we obtain 

(12)' (h' l2>2 h' 21 ) >x + ^2^,2(12)^ = g^ 1 ((12) (hu^i) tl + ^12/121,2(12) 
(21)' (/4 A) i2 + ^12,1 (21); 2 = 92 1 ((21) (h 2l>1 h 12 ) 2 + h 2l h l2>1 (21) 

Therefore we have 

[12]' = g?[l2\ gi (3.18) 

[21]' = g 2 l [21]g 2 (3.19) 

where we define matrix covariants [12] G M mimi and [21] G M mam2 as follows 

[12] = (12) (h 12 ,2h 21 ) A + h 12 h 21j2 (12) A (3.20) 

[21] = (21) (ViM, 2 + ^2i/ii2,i(21) )2 (3.21) 

Before we move to the case n = 3, we compare our covariants (12), (21), [12], [21] with 
Konopelchenko's covariants ( I3.2p : h = a tX +ba—c, k = b^+ab—c, where h and k are covariants 
for the hyperbolic system z xy +az x + bz y +cz = 0. This corresponds to mi = m 2 in the current 

context. As we already know this system can be written in a differential operator form as 

Lz = (d x d y + ad x + bd y + c) z = 0, where the differential operator L = d x d y + ad x + bd y + c 
can be written as 

L = (d x + b)(d v + a)-h 
= (d y + a)(d x + b)-k. 

Therefore, we can rewrite the above system Lz = as 

Au = (3.22) 

Bv = (3.23) 
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where 

For the system (13.221) . we obtain covariant relations: 

(12) = (21)=/i 
[12] = hh xy 
[21] = h x h y 

where mi = m 2 = m and di = d x , d 2 = d y . 
We can easily see that 



[12] + [21] = 1 (hhy) 



Thus, we have 



IV[12]+Tr[21] = iprr(12)(21)]^ 

Similarly, for the system f )3.23p . we have the following relations: 

(12) = (21) = k 
[12] = kyk x 
[21] = kk xy 



These relations give us 



Once again, we have 



[12] + [21] = \ (kk x 



Tr[12]+Tr[21] = i[Tr(12)(21)L 







hi 3 








h 3 i 


h 3 2 


d 3 



The case n = 3: 

Here we consider a differential matrix operator L such that 



L 



where hij G M m<m . = 1, 2, 3) are functions of x\, x 2 and x%. 
Applying the gauge transformation 

L m- L' = g^hg, 

where 

gi(x 2 ,x 3 ) 
#2(^1,^3) 
g 3 (x 1 ,x 2 ) 
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gives us 

Ka = 9T X hi292: h' 13 = gi l h 13 g 3 
Ki = 92 lh 2i9i, K 3 = g2 1 h 2 sg 3 

Ki = g^higi, K2 = 9 3 l h 32 g2 

where g±, g2 and g 3 are invertible square matrices such that g\ G M mimi 

g 3 G M m3m3 . 

By doing some algebraic calculation over the above relations (13.331) - (I3.35p . we obtain 
the following matrix covariants: 



(3.33) 
(3.34) 
(3.35) 

, g 2 G M m2m2 and 



'1 2) - 


'tl2' t 21) \ ± ' } J 


— "'13"'31 




f3 3fi"l 


[23) = 


h 23 h 32 , (21) 


= ^21^12 




(3.37) 


;si) = 


h 31 h 13 , (32) 


= /?-32^23 




(3.38) 


;i23) 


= h 12 h 23 h 31 , ( 


132) = h 13 /i 


32^-21 


(3.39) 


;23i) 


= h 23 h 31 h 12 , ( 


213) = h 21 h 


13^32 


(3.40) 


;3i2) 


= h 31 h 12 h 23 , ( 


321) = « 


21^13 


(3.41) 


12] = 


(12) {h 12 , 2 h 2l 


),i + h 12 h 21>2 


12 ),1 


(3.42) 


13] = 


(13) (h 13>3 h 31 


) A + h l3 h 31>3 


;i3),! 


(3.43) 


21] = 


(21) (h 21A h 12 


) 2 + h 21 h 12>1 


:2i), 2 


(3.44) 


23] = 


(23) (h 23 , 3 h 32 


) )2 + h 23 h 32)3 


:23), 2 


(3.45) 


31] = 


(31) (h 31jl h 13 


) 3 + h 31 h 13A 


:3i). 3 


(3.46) 


32] = 


(32) (/l32,2/i23 


) 3 + h 32 h 23t 2 


:32), 3 


(3.47) 



where (ij), (ijk), [ij] G M miirii (i,j,k G {1,2,3}). 

The question of functional relations between covariants is more subtle than for invariants. 

The general case: 

Let us consider the following differential operator 



L 



/ l mi di 

h>2l Im 2 d2 



\ h n i 



h 



o2 



h2n 

Id) 



where the hij are functions ofx\, x 2 ,--.,x n and the I mi are unit matrices such that h%j G M„ 
and I m . G M minii where i, j G {1, 2, . . . , n}. 
isformation 

L ^ L' = g~ l 1.g, 



The gauge transformation 
where 

/ gi(x2,x 3 , ■ 


V o 



(3.48) 



g 



x n ) 
£2(^1,23, 








\ 



(3.49) 



g n (x 1 ,x 2 , . . • ,x n _i) / 
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gives us 

htj h'ij = g~%jgj (3.50) 

where the gi are square matrices such that G M miTni . 

The relations (I3.50p gives us the following matrix covariants: 

[ij] = (ij)(hij,jh ji ) i + h i:j h ji>j (ij) ji 

(*1*2*3 • • • in) hi 1 i 2 hi 2 i^ ■ ■ ■ ^>i n i\ 

where [ij] e M miJni and (i^is ■■■in) e ^m n m n - 

4. Conclusions and comments 

In this paper, we have dealt with general hyperbolic systems hz = 0. We have used 
a suitable diagonal gauge matrix g, chosen so that it kills diagonal terms ha where i = 
1,2,..., n. We have also obtained the complete set of invariants for general hyperbolic 
systems where rank equals dimension by using the gauge transformation L y L' = g~ Xl Lg. 
Further, we have shown the completeness of a set of simple invariants (reduced invariants). 
We have proved that these invariants form a minimal complete set. 

We have also considered hyperbolic systems hz = where the entries hij are matrices. 
In this case, we are interested in covariants. We have obtained matrix covariants for the 
differential operator L under the gauge transformation. Here we have examined the case 
where rank exceeds dimension. The canonical form of L still suffices where ha = and 
are rectangular matrices. The reduced covariants have been presented but it has not been 
shown that their invariant traces form a complete set. For example, in the case when n = 2, 
we ask the question: Do the covariants (12), (21), [12] and [21] form a complete set? The 
answer depends on the existence of g(xi,X2) so that when 

(12)' = g?(12) gi 
(21)' = g 2 1 (21)g 2 

[12]' = g?[12] gi 

[21]' = g?[21]g 2 

are given then g must satisfy the relation 

g-^g = U. 

The square matrices (ij), (ij)', [ij], [ij] are thus similar (14. II) and so possess as equal 
invariants the traces, say, of their powers: I p = Tr(ij) p etc. But equality of such invariants 
is not sufficient for gauge equivalence of L' and L. There are also invariants associated 
with polynomials in (ij) and [ij] since gi(ij)[ij\ = (ij) [ij]'gi etc, namely, traces of such 
polynomials (cf. (13.26P ). 

Two questions arise for further study: 

(1) What relations on the invariants of these general systems correspond to specialisations 
of L such as self-adjointness? 

(2) Can we establish the existence of a complete, minimal set of trace polynomial invari- 
ants for the systems of Section [3p 
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(3.51) 
(3.52) 
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